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ABSTRACT 

Context. Titanium dioxide, Ti02, is a refractory species that could play a crucial role in the dust-condensation sequence around 
oxygen-rich evolved stars. To date, gas phase Ti02 has been detected only in the complex environment of the red supergiant VY CMa. 
Aims. We aim to constrain the distribution and excitation of Ti02 around VY CMa in order to clarify its role in dust formation. 
Methods. We analyse spectra and channel maps for Ti02 extracted from ALMA science verification data. 

Results. We detect 15 transitions of Ti02, and spatially resolve the emission for the first time. The maps demonstrate a highly clumpy, 
anisotropic outflow in which the Ti02 emission likely traces gas exposed to the stellar radiation field. An accelerating bipolar-like 
structure is found, oriented roughly east-west, of which the blue component runs into and breaks up around a solid continuum 
component. A distinct tail to the south-west is seen for some transitions, consistent with features seen in the optical and near-infrared. 
Conclusions. We find that a significant fraction of Ti02 remains in the gas phase outside the dust-formation zone and suggest that 
this species might play only a minor role in the dust-condensation process around extreme oxygen-rich evolved stars like VY CMa. 
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1. Introduction 


The dust-formation sequence in the outflows of oxygen-rich 
evolved stars is not well understood. It is essential to address 
which gas-phase species provide the primary seeds. Ti 02 is con¬ 
sidered an important seed refractory species with possibly higher 
nucle ation rates than SiO (e.g. |Jeong et aT] |2003[ |Lee et alT} 
2014| ). Moreover, presolar Ti02 grains were tentatively identi- 
fled by [Nittler & Alexand^ ( |1999| ). Since SiO nucleation was 
recently indicated to be more relevant than previously thought 
under the relevant pressure and temperature conditions ( |Nuth &| 
Eerguson[|20Q6t|Gail et al.[[2013| ), it is crucial to characterise the 
role of Ti02. The effect of non-stationarity on the nucleation is 
unknown. Shocks are, however, known to be present in the upper 
atmospheres of these evolved stars (e.g. |Chiavassa et ari|2Qll| ). 


Emission from gas phase Ti02 has only been detected to- 
wardsVY CMa ( [Kaminski et al.[ |2013a|b| hereafter |K-rl3a| 
|K-rl3b| ) . The circumstellar environment of this red supergiant (at 
1.2 kpc; [Zhang et al. | [?0 1 2| ) exhibits a high degree of morpholog¬ 
ical complexity, from optical to radio wavelengths and on spa- 


tial scales from a few to several thousand AU (e.g. [Humphreys 

et al. 

2007 

Kaminski et al.|[2013bl Monnier et al.|[20141 [Muller 

et al. 

r- 

o 

o 

(N 

Shenoy et al.[ 20 13[ [Smith et al.| 20011 Ziurys et al.| 

2007 

). Recent ALMA observations spatially resolve H 2 O maser 


emission, leading to the most accurate determination of the stel¬ 
lar position ( [Richards et al. [2014[ hereafter [R-r 14[ ). [O’ Gorman[ 


[et al.[ ( [2015[ hereafter [Q-rl5[ ) describe the submillimeter contin¬ 


uum emission and report on a bright component south-east of the 
star, indicative of anisotropic mass loss. 

Ti 02 is expected to be consumed by the nucleation process, 
but the detection of Ti02 emission by [K-rl 3 a| already suggested 
that a significant amount could possibly survive the dust forma¬ 
tion. [K-rl3a[ did not spatially resolve Ti02 emission at angular 
resolutions >1". The ALMA observations now allow us to char¬ 
acterise the emission and the role of Ti 02 in the dust-formation 
process in more detail. 


2, Observations 


We retrieved ALMA science verification data on VY CMa from 
the ALMA archives. The o bservati ons and data calibration and 
reduction are described by [R-fl4[ Table shows the spectral 
coverage and representative rms noise values for the six spectral 
windows in ALMA’s band 7 (~0.9 mm; ~320 GHz) and one in 
band 9 (-0.45 mm; -660 GHz). With projected baselines of 14 m 
up to 2.7 km, the spatial resolution at -320 GHz and 658 GHz 
is - 072 and - 071, respectively, and the maximum recoverable 
scales are 873 and 470. The synthesised beam sizes are obtained 
using natural weighting. 

The data reduction and quality of the continuum emission 
ima ges at 3 21 GHz and 658 GHz are discussed in detail by [R-fl4 
and |0-fl5[ No imaging artefacts are expected to arise from the 


array configuration owing to the excellent coverage of the visi¬ 
bility plane. 
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Table 1. Spectral coverage of the ALMA observations. 


Setting 

Frequency range 
(GHz) 

Resolution 

(MHz) 

Rms noise 
(mJy beam“^) 

Beam 

("x",P.A.) 

321 

309.501 -310.437 

0.98 

3 

0.23 x0.13,28.27° 

321 

310.501 -311.437 

0.98 

2 

0.23 x0.13,28.48° 

325 

312.035 - 313.908 

1.9 

3 

0.24x0.13,29.48° 

321 

320.725 -321.662 

0.98 

3 

0.22x0.13,27.31° 

321 

321.986-322.923 

0.98 

3 

0.22x0.13,27.55° 

325 

324.181 -326.054 

1.9 

6 

0.22x0.12,29.27° 

658 

657.002 - 658.873 

3.9 

15 

0.11 X 0.06,29.89° 


Notes. Setting names refer to those used for the released ALMA CSV data, focussing on the H 2 O masers at 321 GHz, 325 GHz, and 658 GHz. 


We note that the spectral window 312 - 314 GHz suffered 
from poor continuum subtraction owing to line crowding and a 
consequent lack of line-free spectral range. This causes a frac¬ 
tion of the continuum emission to leak into the channel maps of 
the Ti 02 lines at ~312GHz (Figs. [A^ |A.5| [Al 6 | ). The emission 
seen at the position of clump C can be entirely attributed to the 
continuum itself, i.e. we can rule out molecular contribution at 
this position. The lines presented in Fig.l^and Table|^ and con¬ 
sequently also the reported peak and integrated intensities, were 
corrected for this effect on a line-to-line basis. 

The imaged Ti 02 lines in the 312-314 GHz spectral window 
also show a contribution north-north-east and south-south-west 
of the star. These features are artefacts of the cleaning procedure 
likely caused by imperfect phase corrections on a number of 
intermediate baselines. The phase corrections were transferred 
from the self-calibrated 325 GHz maser line located in the atmo¬ 
spheric absorption region. 


3. Results 


We analysed spectra extracted for a 1" diameter region around 
the stellar position. We show below that no Ti 02 emission is de¬ 
tected beyond this aperture. Ti 02 identifications are base d on the 
Cologne Database for Molecular Spectroscopy (CDMS; [Muller 


jet al.||200T]|2005t|Brunken er3Il|2008| ). We detect 15 lines with 
upper-level energies E^^/k in the range 48 - 676 K and signal-to- 
noise ratios S/N ^ 5-11 aX velocity resolutions 0.9 - 7.6 km s 


(Table 13 Fig. [ 3 . Of the Ti02 lines detected by |K-rl3a||K+m 


only those at 310.55 GHz and 310.78 GHz are observed with 
ALMA. The peak fluxes of the ALMA and SMA detections 
are consistent within the uncertainties. Moreover, of the other 
13 lines of Ti02 detected with ALMA, none were detected in 
the SMA survey, owing to a noise level in the ALMA data that 
was approximately 10 times lower. We detect no Ti02 emission 
at ~ 660 GHz owing to the higher rms noise (Table [ 3 . We de¬ 
tect no isotopic variants of Ti, consistent with the solar isotopic 
abundance ratios. 

We expect no large flux losses in the ALMA observations of 
Ti02, given the ~8" recoverable scale at -320 GHz. The pres¬ 
ence of large-scale emission with low surface brightness cannot 
be entirely excluded, but is unlikely and will not change our main 
conclusions. 


Blending, proximity to strong lines, line crowding, and in¬ 
trinsic line-shape irregularities complicate the identification of 
Ti02 emission. The lines at 310.55 GHz and 311.46 GHz are 
only partially covered in the observations; we believe the former 
to be a firm detection and the latter to be tentative. In both cases, 
no other candidate could be identified. Uncertainties on the rel- 
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evant Ti02 frequencies are < 2 MHz, except for the transition at 
311.462 GHz, where it is 5.555 MHz ( [Brunken et ari[2008l ). 

The Ti02 lines detected with ALM A exhibit broad line pro¬ 
files, as do those detected by K-rl3a Assuming a stellar ^lsr 
of 22 km s“^ we find that the emission is very asymmetric in ve¬ 
locity space, ranging typically from —15 km s“^ to ~ 60 km s“^ 
(Fig. [ 3 . The apparent extension of the blue wing of some lines to 
—45 kms“^ is likely coincidental and due to line blending, but 
we cannot exclude an actual TiO? contribution. Blend candidates 
are indicated in Fig.[^and Tablej^ The lines with E^^/k = 48 kR 
58 K, 182 K, 532 K exhibit features in their red wings that extern 
to fLSR - 105 km s“^ A similar high-velocity outflow is also vis¬ 
ible for some lines presented by |K-rl3^ 


3.1. Spatial distribution 

The detected Ti02 emission is, for the most part, spatially re¬ 
solved by the ALMA observations with a 0723 x 0713 beam 
(see Figs. |A.lf|A.16| ), but several transitions show spatially un¬ 
resolved emission at a 2kms“^ velocity resolution. We note 
a very complex behaviour of the emission peaks, with some 
only appearing in one channel and not in the neighbouring 
channels. The often spatially unresolved (<145 AU) peaks in 
the different velocity channels of all lines imply a clumpy 
and/or anisotropic wind. The maps for 310.55GHz (£’up/^ = 
181 K) and 310.78 GHz {E^^jk = 182K) reveal that these two 
lines behave almost identically. Corresponding velocity chan¬ 
nels (Figs. [ATT||A.2| ) trace roughly the same regions of the cir- 
cumstellar environment and exhibit similar intensities, imply¬ 
ing that these two lines are very tightly coupled in their exci¬ 
tation. Figure [3 shows the morphology of the Ti02 emission 
at 310.78 GHz and the positions of the star (VY), the contin¬ 
uum component C ( |Q-rl5[|R-rl4| ), and the so uth-west clump de- 
tected at 4 ~ 1 - 5 pm at ~ 1" from the star ([Smith et n^|2001t 


[Humphreys et al.[[2007t[Shenoy et al.[[2013|). 

Given the similar energy levels, quantum numbers, and 
Einstein-A coefficients of the transitions at 310.55 GHz, 
310.78 GHz, and 324.49 GHz, one expects similar line intensi¬ 
ties and spatial distributions. However, although th e chan nel-to- 
channel peaks of the 324.49 GHz emission (Fig. |A.11[ ) corre¬ 
spond quite well to those of the other two, which strengthens its 
identification as Ti02, its intensity is clearly lower. It is hard to 
explain this discrepancy, but we note that as a consequence of 
the sensitivity of the atmospheric transmission at ~ 325 GHz to 


See 


K+13b 


for a discussion of Ulsr- 

subtraction is hampered by the high line density. The red 


^ Continuum s 

wing of the line is affected, but still shows the feature they have in com¬ 
mon. 
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Fig. 1. Ti02 spectra extracted for a 1" diameter aperture around the stellar position. The vertical dashed lines indicate the stellar tiLSR of 22 km s ^ 
the shaded areas the [Ilsr- ranges from TableWe indicate identifications of species other than Ti02 in the panels. 

Table 2. Overview of detected Ti02 lines. Parameters listed for spectrum extracted in a 1" diameter aperture centred on the star. See also Fig.[^ 


Transition 

Uab 

(MHz) 

Eaplk 

(K) 

fmin t’jYiax 

(kms“^) 

Ai; 

Rms 

(mJy) 

Peak 

(Jy) 

/ 

(Jykms“^) 

Size 

(") 

sw 

tail 

22(1,21)-21(2,20) 

310554.735 

180.5 

-15.2 

51.3 

1.0 

16 

0.22 

9.57 

0.9 

Y 

23(1,23)-22(0,22) 

310782.713 

182.4 

-43.9 

102.3 

1.0 

16 

0.27 

15.28 

0.9 

Y 

40(8,32)-39(9,31) 

311462.082 

675.9 

19.4 

66.7 

1.9 

10 

0.10 

2.31 

0.8 

- 

7(5,3) - 6(4,2) 

312248.341 

47.9 

-43.7 

80.96 

1.9 

16 

0.29 

20.32 

0.8 

N 

10(4,6) - 9(3,7) 

312732.066 

57.7 

-19.4 

105.1 

1.9 

16 

0.15 

9.43 

0.5 

? 

30(3,27)-30(2,28) 

312816.809 

354.4 

-68.2 

60.0 

1.9 

16 

0.16 

12.98 

0.5 

? 

11(4, 8)- 10(3,7) 

321401.936 

65.7 

-33.9 

50.5 

0.9 

24 

0.19 

7.24 

0.8 

? 

28(2,26)-28(1,27) 

321501.043 

298.8 

-35.0 

51.2 

0.9 

24 

0.13 

4.76 

0.9 

N 

35(8,28)- 35(7,29) 

322333.594 

532.4 

-36.4 

108.1 

1.0 

19 

0.15 

7.03 

0.8 

N 

33(8,26)- 33(7,27) 

322612.696 

481.4 

-0.1 

63.9 

1.0 

18 

0.16 

4.18 

0.9 

N 

23(2,22)-22(1,21) 

324492.930 

196.1 

-47.9 

79.4 

1.9 

29 

0.21 

9.64 

<0.5 

N 

37(8,30)- 37(7,31) 

324965.693 

586.5 

3.6 

59.8 

1.9 

30 

0.15 

4.92 

0.5 

N 

26(1,25)-26(0,26) 

325322.035 

246.5 

-25.2 

81.8 

7.6 

21 

0.11 

1.85 

0.5 

N 

26(4,22)-25(5,21) 

325500.415 

281.5 

-7.7 

55.7 

3.8 

28 

0.13 

2.75 

0.5 

N 

28(8,20)-28(7,21) 

325600.792 

367.0 

-0.4 

90.2 

3.8 

28 

0.18 

3.89 

0.5 

N 


Notes. Columns are transition, rest frequency, upper-level energy, minimum and maximum jilsr reached, velocity resolution and rms noise at 
which the identification was made, peak and integrated fiux, maximum diametric size of the integrated emission at 3cr level, and a marker for 
emission in the south-west tail (Y, N, ?, and - for yes, no, maybe, and unknown). Species other than Ti02 detected inside the plotted spectral 
windows are marked in Fig.[2 


the atmospheric water-vapour content the rms noise is ~ 4 times 
higher and the flux calibration could be compromised. 

The emission at 322.61 GHz, 322.33 GHz, and 324.96 GHz 
behaves similarly but with weaker and spatially somewhat more 
conflned red-wing emission. Unfortunately, many of the Ti02 
lines are blended (both in frequency and spatially) with emission 
from other species. 

The discussion below is focussed on the morphology at 
310.78 GHz with Fig. as a visual guide. Overall, the detected 
Ti 02 emission moves from west to east, across the stellar po¬ 
sition, with velocities evolving from reddest to bluest. All lines 
show emission close to the star at the stellar ulsr- In particular, 
all lines with F’up/^ > 360 K are centred on VY and emit mainly 
across l>lsr ~ 0 - 40kms“^ (Fig.[^, implying that these trace a 
central part of the outflow that may be accelerating. 

At l>lsr ^ 44kms“^ the emission is mainly situated west of 
VY with a marginal contribution close to it. The emission is 


highly variable with ^lsr, doubly peaked in most velocity chan¬ 
nels, and shows a hook-like feature at its western e dge. The lat- 
ter could relate to the north-west knot deflned by [Humphrey^ 
et al.| ( |2QQ7l ) based on HST observations, though more analysis 


is needed to investigate the nature of this feature. 

At 22 < ^lsr < 44kms“^ we And multiple emission peaks 
in most velocity channels, with those closest to VY slightly 
brighter than the west offset peaks. The integrated emission is 
elongated along a direction roughly parallel to the axis connect¬ 
ing the peak position of C and VY, at a RA. 125°. 

At 18 < t>LSR ^ 22kms“^ the transitions at ~310GHz show 
a bright tail extending south-west of VY at a RA. 220°, i.e. 
almost perpendicular to the axis connecting VY and C. This tail 
reaches ~ 1" away from VY, out to the south-west clump of 


IShenoy et al. 

(|2013) and agrees very well with the features de- 

tected by e.g. 

Smith et al. (|2001|) at wavelengths T ~ 1 -2.14 jum. 


The emission at 310.78 GHz extends slightly further south-west 
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Fig. 2. Ti02 morphology. Colour maps of emission at 310.78 GHz integrated over the t^LSR-ranges indicated at the top left of each panel, cut off at 
3cr. Red contours s how the 321 GHz continuum at [3,20,40,60,80]cr; green contours show HST emission at [3,5,7,10,20,30,40,50,100, 200]cr 
( Smith et ari|2001| ). In the first panel we mark the position of the star (+, VY; black) and of the continuum component ( x, C; red) to the south -east 
( 0+15||R+14| ), and the position and approximate extent of the south-west clump (SW, dashed 1" diameter circle; red |Shenoy et al.|[2013| l. The 
apparent north-south emission is thought to arise from dynamic-range limitations in the peak channels. 


than that at 310.55 GHz. Since the most extended emission ex¬ 
ceeds 5(T, we suggest that the difference is real and that dumpi¬ 
ness in the outflow strongly influences the excitation of individ¬ 
ual lines. 


At 312.73 GHz and 312.82 GHz we And signs of emission in 
the south-west tail, but these are likely artefacts from the clean¬ 
ing procedure (see Sect. |^. At 311.46 GHz we see emission 
within the spatial region where the wouth-west tail is located, 
although at higher ^lsr than the transitions at -310 GHz. This 
could be due to the south-west tail covering a wider velocity 
range than reported above, by a line blend, or a misidentiflcation 
of the line as Ti02. We unfortunately do not have information 
for l>lsr < 22 km s“^ for this line. No other transition shows de¬ 
tectable emission within the south-west tail, but most show a 
slight bulge of emission at - 0'.'2 south-west of the star, at the 
base of the south-west tail. Figure [ATTT] shows a comparison for 
all lines to the south-west tail observed at 310.78 GHz. 

AtfLSR < 18kms“^ the high-intensity emission is situated 
entirely east of VY, elongated, and oriented at - 15° - 25° east 
from north. With bluer velocities the emission moves towards 
C and then appears to break up with a northern peak brighter 
than the southern one (e.g. Fig. |A.2| ). Remarkable is that the 
low-intensity component of the 310.78 GHz transition closely 
resembles the scattered light at 1 pm. We And this strong cor¬ 
respondence at these blue velocities for no other Ti 02 transition. 
If the high- and low-intensity components have different intrin¬ 
sic wind velocities, they could be spatially separated and trace 
different parts of the outflow. 


13 

4 
d 

d 11 

o 

10 

0 200 400 600 800 

Eu,/k [K] 

Fig. 3. Rotational diagram. Lines indicated in red are blended or only 
partially covered in the observations and are excluded from the fitting 
procedure. All intensities, source sizes, and rms noise values are taken 
from Table [2] The fit results and uncertainties are shown with the blue 
and grey dashed lines and are indicated at the top right. 


= 


V = 5.65+/- 1.33 10^= cm^ E 
198.0 +/- 28.5 K 




Th e high dipole moment of Ti02 (6.33 Debye; [Wang et al.[ 
2009|) supports efficient radiative excitation and T^ot could hence 


reflect an average continuum brightness temperature as opposed 
to a gas kinetic temperature in the case of collisional excitation. 
Additionally, the large dipole moment induces electron-Ti02 
collision rates large enough to exceed H 2 -Ti 02 collision rates 
if the fractional ionisation exceeds a few 10“^. We could not And 
a value for VY CMa, but ba sed on the result ^e/ ^H = 3.8 x 10“"^ 
for the red supergiant a Ori ( [Harper et al.[|200l] ) this type of ex¬ 
citation could be relevanj^ We discuss the competition between 
collisional excitation and radiative excitation in more detail be¬ 
low, in Sect. [4.1 j 


3.2. Excitation conditions 


4. Discussion 

4.1. Excitation of titanium dioxide 


We investigate the excitation of Ti02 lines via a rotational dia¬ 
gram analysis (Fig. [^. Intensities, source sizes, and rms noise 
values are taken from Table Severely blended or only par¬ 
tially covered fines are excluded from the analysis. We derive a 
source-averaged column density Acoi = 5.65 + 1.33 x 10^^ cm“^ 
and a rota tional te mperature Yj-ot = 198.0 + 28.5 K, in agree¬ 
ment with |K-rl3a[ We note that the kinetic temperature in the 
excitation region of Ti02 varies from more than 1000 K down 


to -100 K (Decin et al.l 

. 2006 

). Assuming an avera 

ge mass-loss 

rate of 2 x 10“"^ M© yr“ 

' (e.g. 

De Beck et al. 2010 

), an average 


find an abundance Ti02/H2 3.8 + 0.9 x 10 


Many effects need to be taken into account in a full treatment 
of the excitation and radiative transfer of rotational transitions 
of Ti02. Important molecular data are lacking: although the ro¬ 
tational energy levels and radiative data are well determined for 
Ti02 in CDMS ( [Muller et al.|[200T][20051 [Brunken et al.|[2008l ), 
no rates exist for excitation of Ti 02 by hydrogen-impact nor have 
the vibration-rotation spectra been fully analysed. 

For illustration of the transition rates for Ti02, consider 
two of the observed rotational transitions at 312.248 and 


VY CMa is of spectral type M2.5-M 5e la ([Houk & Smith-Moore 


1988[ ), a Ori of spectral type Ml-2Ia-ab { [Keenan & McNeil 1989[ l. 
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Table 3. Summary of properties of selected transitions of Ti02. 



y 

Au/ 

Au,{ 

no 

Trad 

Prad 

or vibrational band 

(GHz) 

(s->) 

(s->) 

(cm“^) 

(K) 

(s->) 

7(5,3)-6(4,2) 

312.248 

4.34 X 10-3 

5.04 X 10-3 

5 X 102 

26.2 

5.7 X 10-3 

7(5,3)-7(4,4) 

206.094 

6.95 X lO--^ 



15.3 

7.6 X 10-4 

7(5,3)-8(4,4) 

84.567 

7.31 X 10“® 



7.15 

9.6 X 10-® 

37(8,30)-37(7,31) 

324.966 

4.76 X 10-3 

4.63 X 10-2 

5 X 10^ 

27.7 

6.3 X 10-3 

37(8,30)-36(9,27) 

205.259 

2.62 X 10-‘3 



15.2 

2.9 X 10-4 

37(8,30)-36(7,29) 

821.811 

2.96 X 10-2 



94.3 

5.7 X 10-2 

37(8,30)-36(5,31) 

1235.577 

1.17 X 10-2 



132.4 

2.0 X 10-2 

vi = 1 - 0 

28840 

3.6 

3.6 

a 

350 

7.2 X 10-2 

y2 = 1 - 0 

9863 

7.2 X 10-2 

7.2 X 10-2 

a 

300 

2.0 X 10-2 

U) 

II 

1 

0 

28030 

42 

42 

a 

350 

9.0 X 10-' 


Notes. Columns list the transitions (rotational and vibrational), the frequency v at which these transitions occur, their transition probability A^/, 
inverse lifetime critical density hq, radiative temperature Trad and the calculated pumping rate prad- Appropriate collision rates are 

unknown; therefore, no critical density is tabulated for the vibrational transitions. 


324.965 GHz with properties summarised in Table[^ Collisional 
excitation at the kinetic temperature of the gas can dominate only 
when the downward rate of collision-induced transitions greatly 
exceeds the downward rates of radiative transitions for each up¬ 
per state. This condition can be translated to a critical density 
^0 = where is the spontaneous transition prob¬ 

ability of each transition from upper state u to lower state To 
estimate this density, we assume a characteristicj^ quenching rate 
coefficient qo = 10“^^ cm^ s“^ for rotational transitions induced 
by collisions with neutral species H or H2. Values for no are listed 
in Table[^ Assuming a mass-loss rate of 10“"^ M© yr“\ such den¬ 
sities are however reached only at distances from the star lower 
than ~0''065, whereas Ti02 is excited over a much larger region, 
where the density quickly decreases to a few 10^ cm“^. 

In addition, collision rates for electron impact have been 
computed in the Born approximation. Because the electric dipole 
moment of Ti02 is so large, electron-impact rates are expected 
to show a strong propensity for radiatively allowed transitions. 
Computed rates are likely to be accurate within 50%. Under 
an assumed very low fractional electron density of 10“^, the 
collision-induced downward transition rate for the 312 GHz tran¬ 
sition is for 5% due to electron collisions. Assuming the frac¬ 
tional ionisation reported for aOri (3.8 x 10“"^) electron colli¬ 
sions would completely dominate over neutral (hydrogen) colli¬ 
sions by a factor of the order of 100 or more and might compete 
well with infrared pumping. 

We note that the upper state of the 312 GHz transition decays 
mainly by the observed transition itself, while the upper state 
of the high-excitation 324 GHz transition is depopulated more 
rapidly by submm-wave transitions at 821 and 1235 GHz. As a 
consequence, a much higher density would be required to excite 
the 324 GHz line by collisions than the 312 GHz line. 

The continuum intensity of VY CMa is so large at in¬ 
frared and submm wavelengths that radiative excitation (pump¬ 
ing) must be taken into account. If the observed continuum flux 
of VY CMa is assumed to arise within a 0'.'3 diameter region 


Such a value is typical of the largest downward collisional rates 
for H 2 collisions with a heavy, polar mole cule. For example , accu- 
rate collision rates have been computed by [Cernicharo et al.| ( |2011| ) 
for H 2 on SO 2 , a heavy molecule with a relatively large dipole mo¬ 
ment (1.63 Debye). The SO 2 quenching rates at low temperature are 
~ 2 X 10-1^ 


( |Q-fl5| ) and to be diluted by a geometrical factor 1/9 over the 
0'.'9 extent of the observed Ti02 emission, then the average 
Planckian radiation brightness temperatures are estimated to be 
Trad ^ 25 K at 312 to 324 GHz. The corresponding pumping 
rates (absorption and stimulated emission) in this radiation fleld 
can be expressed as 

_ 

exp(/zy/^rrad) - 1 


for each radiative transition that connects states of interest. The 
basic radiative data for Ti02 are collected in Table Although 
the fundamental bands of the three vibrational modes have not 
been fully analysed rotationally, the band frequencies and band 
strengths are approximately known ( |Grein[|2007) ). In the adopted 
continuum model, the infrared intensity is high enough to drive 
absorption in vibration-rotation lines at rates of the order of 0.01 
to 1.0 s“^. In order for collisional excitation in a pure rotational 
transition u ^ iio compete with infrared pumping in VY CMa, 
a density much greater than 


«o ~ (p(vi) + P(V2) + P(V 3 ) + ^ Au,t)lqQ « 10'° cm ^ 

( 

would be required. The radiative rates in Table [^suggest that ra¬ 
diative excitation is likely to be very important for Ti02 in VY 
CMa. Therefore, the rotational temperature derived in Sect. |3.2| 
might have no direct relationship to the kinetic temperature. Be¬ 
cause the observed transitions in the present study span a wide 
range of excitation energies, the simple rotation-diagram analy¬ 
sis still provides a useful first estimate of the molecular column 
density and abundance. 


4.2. Outflow components 


From the comparison to the continuum at ~ 321 GHz, we derive 
that Ti02 is excited in the directions with lower dust densities. 
The absence of detected Ti02 emission north of the star could 
then point to efficient obscuration of this part of the outflow, in 
line with the observations of e.g. [Smith et al.| ( [200T] ). Attenuation 


of the stellar radiation fleld to the east and west of the star is 
limited. From this, we do not expect an equatorial enhancement 
of the mass-loss rate, since this would likely have induced more 
efficient dust formation, which is not seen in the continuum. 
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The Ti02 emission traces multiple wind components. We 
find a red outflow to the west and a blue outflow to the east. 
With the clear exception of the interaction of the Ti02 gas with 
clump C in the east, the two seem roughly symmetric around the 
star and aligned with the axis connecting VY and C. We rule out 
an equatorially enhanced environment such as an expanding disk 
or ring, based on the spatial distribution of the Ti02 emission at 
different ^lsr- We rather suggest an accelerating bipolar-like out¬ 
flow at lower densities. We also find a predominantly blue south¬ 
west outflow, connecting the star and the south-west clump, ap¬ 
proximately perpendicular to the VY-C connecting axis. We 
find no north-east counterpart in Ti02 emission, likely imply¬ 
ing that the south-west tail is indeed caused by an event in one 
preferred direction, as opposed to a bipolar event. 

4.3. Interaction of blue outflow with clump C 

Whereas the H2O maser emission in the “valley” between VY 
and C implies that C is close to or in the plane of the sky ( |R-rl4| ), 
the observations suggest that the Ti02 gas breaks up around C 
while moving towards the observer, placing C - at least partially 
- in front of the plane of the sky. We therefore deem it likely 
that the H2O masers and the Ti02 emission probe parts of the 
outflow east of VY with different physical properties. Whereas 
the masers are probably excited through shocks at high densities, 
Ti02 is more likely excited through radiation, at lower outflow 
densities. In the denser regions, Ti02 might not be excited and/or 
it might be efficiently depleted from the gas phase. The latter 
is, however, less likely (see below). We therefore suggest that 
Ti02 traces the blue-shifted wind to the east of VY with lower 
densities which runs into and curves around C. 


4.4. Titanium dioxide and scattered light 


K-rl3b| reported emission offset by ~ 1" south-west from the 
central molecular emission for multiple species. However, ow¬ 
ing to lower sensitivity they did not find indications for this in 
the Ti02 emission, whereas we clearly detect the south-west tail 
at-310 GHz (see Sect.[3T]). 

The likelihood of radiative excitation of Ti02 and the agree¬ 
ment between the south-west tail in the Ti02 emission and the 
scattered stellar light at 1 pm from [Smith et ah] ( |2001| ) suggest 
the presence of a tail of gas and dust lit up by a stellar radia¬ 
tion field less attenuated than at other position angles. The bulge 
south-west of the star seen for most Ti02 transitions is consis¬ 
tent with this. We infer a lower degree of attenuation and from 
that a lower (dust) density south-west of the star. We speculate 
that the suggested localised ejection leading to the formation of 
the south-west clump ( [Shenoy et al.||2013] ) could have created a 
cavity in the dense material close to the star. 

Given the agreement with the blueshifted Ti02 emission, the 
outflow traced by the scattered light is likely oriented out of 
the plane of the sky. The lack of consistent ^lsr coverage over 
the different Ti02 transitions, however, further complicates the 
three-dimensional and kinematic constraints. In-depth analysis 
of all other emission lines detected with ALMA is needed to 
constrain the basic properties of the south-west tail. 


4.5. Titanium dioxide and dust formation 

If the nucleation of Ti02 can occur at < 2000 K (e.g. ILee et al. 


close to the star, i.e. within a few stellar radii. From the ALMA 
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2014), one expects Ti02 depletion from the gas phase initiatec 


observations |Q-r 15| put an upper limit to the dust-condensation 
radius of ~ 0'.'06 + 0'.'02 10 stellar radii. The Ti02 emis¬ 

sion is present close to the star and out to <0'.'45 in directions 
without extended continuum emission. Where Ti02 and the bulk 
dust coexist, the emission extends to ~0'.'20, well beyond the 
dust-condensation radius. Since the Ti02 emission is seen at ra¬ 
dial distances from the star far beyond the dust-condensation re¬ 
gion, and along directions where dust continuum is observed, we 
claim that Ti02 is not a tracer for low grain-formation efficiency. 
The strong correspondence between the Ti02 emission and the 
south-west tail of dust-scattered stellar light (see Sect |3.1| and 
Fig. 1^ supports this. We also considered the high derived col¬ 
umn density and suggest that Ti02 plays only a minor role as a 
primary dust seed around VY CMa. 

Since the ALMA observations covered no transitions of TiO, 
we cannot expand the discussion of |K-rl3a] on the relation be¬ 
tween TiO, Ti02, and the dust. 


5. Conclusion 

We detect 15 transitions of Ti02 in high-resolution ALMA ob¬ 
servations of the red supergiant VY CMa. The main emission re¬ 
gion spans ~0'.'9 (-1080 AU) in a roughly east-west oriented di¬ 
rection, centred on the star. The behaviour of the gas is complex 
throughout t>LSR-space with bright peaks appearing only in very 
confined l>lsr -ranges, implying that the outflow is very clumpy, 
on spatial scales not resolved by the current observations with a 
0'.'23 X 0'.'13 beam. We observe a tail of Ti02 emission extend¬ 
ing out to -1" south-west of the star, consistent with structures 
seen in the optical and near-infrared. It is oriented out of the 
plane of the sky and mainly covers projected velocities of a few 
kms“\ but reaches up to -40kms“^ in some cases. We suggest 
that the Ti02 in this tail is illuminated by stellar radiation pene¬ 
trating through a low-density cavity in the south-west part of the 
circumstellar environment of the star, potentially created by the 
ejection that led to the south-west clump. Within a bipolar-like 
Ti02 outflow, the blue-shifted emission exhibits a strongly dif¬ 
ferent orientation and behaviour than the red-shifted emission, 
suggesting that the stellar wind runs into the large body of dust 
situated -0'.'335 (-400 AU) south-east of the star. 

We suggest that Ti02 might play only a minor role in 
the dust-condensation process in the complex outflow around 
VY CMa, and potentially also around other oxygen-rich evolved 
stars with extreme mass outflows. High-resolution imaging is 
however still needed to correlate the emission of Ti02 with that 
of TiO, and to further investigate the relative importance of sili¬ 
con, titanium, and other metals in the dust condensation. 
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Appendix A: Maps of titanium dioxide emission 


Figures |A.1| to |A.15 
emission lines (Tab! 


show channel maps of the detected Ti02 
e Fig. [T]) at a velocity resolution of 
2kms“\ covering the range -14 < l>lsr ^ 78kms“^. This l>lsr- 
range covers the bulk of all Ti02 emission; emission at more 
extreme velocities is no longer visible in the channel maps. Fig- 
ure |A.16| shows integrated-intensity maps for all listed lines, cov¬ 
ering the t>LSR-ranges indicated in Table Figure [ATTT] shows a 
comparison of the Ti02 line emission to the south-west tail de¬ 
tected at 310.78 GHz, which is discussed in Sects. |3.l|and|4.4| 
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Fig. A.l. Channel maps of the Ti02 emission at 310.55 GHz, at a 2 km s"^ velocity resolution. Black contours show the continuum measured with 
ALMA at 321 GHz { |0+15|[R+14| ). The stellar position is indicated with a white cross. Spatial scales are indicated in the top left panel and are the 
same for all panels. The colour scale starts at the 3cr level and is plotted as the square root of the flux for increased contrast. 
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Fig. A.2. Same as Fig. A.l but for Ti02 emission at 310.78 GHz. 
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Fig. A.3. Same as Fig. A.l but for Ti02 emission at 311.46 GHz. 
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Fig. A.4. Same as Fig. A.l but for Ti02 emission at 312.25 GHz. 
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Fig. A.5. Same as Fig. A.l but for Ti02 emission at 312.73 GHz. 
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Fig. A.6. Same as Fig. A.l but for Ti02 emission at 312.82 GHz. 
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Fig. A.7. Same as Fig. |A.l| but for Ti02 emission at 321.40 GHz. Artefacts in channels with i>lsr > 55kms ^ are due to the presence of strong 
SO 2 emission. See also Fig.[2 
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Fig. A.8. Same as Fig. A.l but for Ti02 emission at 321.50 GHz. 
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Fig. A.9. Same as Fig. A.l but for Ti02 emission at 322.33 GHz. 
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Fig. A.IO. Same as Fig. A.l but for Ti02 emission at 322.61 GHz. 
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Fig. A.ll. Same as Fig. A.l but for Ti02 emission at 324.49 GHz. 
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Fig. A.12. Same as Fig. A.l but for Ti02 emission at 324.96 GHz. 
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Fig. A.14. Same as Fig. A.l but for Ti02 emission at 325.50 GHz. 
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Fig. A.16. Integrated-intensity maps of Ti02 emission. Colour maps show the intensity integrated over the ulsr- ranges indicated in Table and 
Fig. and cut off at 3cr. Contours show the ALM A 321 GHz continuum. Labels indicate the positions of the star (+, VY) and the contin uum 
component (x, C) to the south-east ( |0-hl5[[R-hl4| ), and the position and approximate extent of the south-west clump of |Shenoy et al.| ( [2013| SW 
and a dashed 1" diameter circle). The presence of cleaning artefacts, mainly at ~312 GHz, is addressed in Appendix [A| 



























































































Dec. offset ["] Dec. offset ["] Dec. offset ["] Dec. offset ["] 


A&A-25990_DeBeck, Online Material p 25 


Flux [Jy beam ’ km s '] 


0.005 0.015 0.025 0.034 0.044 0.054 0.063 0.073 



1.5 1.0 0.5 0.0 -0.5 -1.0 -1.5 

R.A. offset ["] 


(a) Ti02 at 310.55 GHz 


Flux [Jy beam ' km s’] 



(b) Ti02 at 312.25 GHz 


Flux [Jy beam ’ km s ’] Flux [Jy beam ’ km s ’] 



Flux [Jy beam ’ km s ’] 



(e) Ti02 at 321.40 GHz. 


Flux [Jy beam ’ km s’] 



(f) Ti02 at 321.50 GHz 


Flux [Jy beam ’ km s ’] 



(g) Ti02 at 322.33 GHz 


Flux [Jy beam ’ km s ’] 



(h) Ti02 at 322.61 GHz 


Flux [Jy beam ’ km s ’] 



Flux [Jy beam ’ km s’] 



(i) Ti02 at 324.49 GHz Q) Ti02 at 324.97 GHz 


Flux [Jy beam ’ km s ’] 



(k) Ti02 at 325.32 GHz 


Flux [Jy beam ’ km s ’] 



Flux [Jy beam ’ km s ’] 



(m) Ti02 at 325.60 GHz 

Fig. A.17. Comparison of Ti02 emission lines listed in Table|^to south-west tail at 310.78 GHz (red contours) and HST image (black contours). 
Emission integrated over 19 < t^LSR ^ 22kms"^ plotted at > 3cr (colour scale). Emission of 311.46 GHz was omitted since the relevant i^LSR-range 
is not covered by the observations. 




























































































































































